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Abstract. I discuss recent experimental data on the Cronin effect in deuteron- 
I gold collisions at the top RHIC energy, in a pseudorapidity range 77 £ [—2,3]. 

^2 ' Two theoretical approaches are compared and contrasted: the pQCD-based 

I Glauber-Eikonal model and Colour Glass Condensate models. Neither can 

ly-j ■ describe the Cronin effect over the whole pseudorapidity interval up to now 

I explored experimentally, its most mysterious and intriguing part being at neg- 

■ ative rapidity. 

P_ 
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^ ' In hadron-nucleus {p + A) and nucleus-nucleus {A + A) collisions at relativistic 

P3 . energy the hadron transverse momentum spectra at moderate pT^2-6 GeV are 

enhanced relative to linear extrapolation irom p + p reactions. This "Cronin effect" 
has been observed on an energy range ^/s w 20 — 200 GeV in both p + A and A + A 
' collisions [1]. It is generally attributed to multiple scattering of projectile partons 

propagating through the target nucleus [2] , and is a sensitive probe of initial-state 
modifications of the nuclear wave function [p + A collisions) and of final state 
in- medium effects {A + A collision) [3]. The Cronin effect may be quantified by 
taking the hadron px spectrum in p + A collision in a given centrality class (cel.), 
normalizing it to binary scaled p + p collisions by the inverse thickness function T^, 
and finally dividing it by the p + p spectrum: 



1 dNP^:^''^, / daPP^''^ 

[cel.] / . 

Ta(c.cL) dq^dy / dq^dy 



To cancel systematic errors as much as possible, it is also customary to take the 
ratio Rep of a given centrality class to the most peripheral one. 
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1. Recent experimental results at RHIC 

The BRAHMS results [4] on charged hadron R^p in d + Au colhsions at ^/s — 200 
GeV A in a rapidity range 77 = — 3.2 are shown in Fig. 1, top panel. Two features 
are apparent. First, the Cronin peak at pT — 3 GeV in the mid-rapidity bin is 
progressively suppressed. At 77 > 2, Rep is always smaller than 1, tending to a 
plateau at larger transverse momenta. Second, the centrality dependence of the 
effect is reversed going from mid- to forward rapidity. 

Peripheral collisions are expected to behave similarly to p + p collisions: the 
more peripheral the collision, the thinner the nucleus, the smaller the rescattering 
probability. Therefore, the Rep ratio should behave similarly to RdAu- Quite strik- 
ingly, this is not the case, as shown in Fig. 1, bottom panel. The peak in RdAu 
is suppressed in going from 77 = to 77 = 1, as confirmed also by PHOBOS data 
[5], and looks similar to Rep. However, at 77 > 2 the suppression doesn't seem 
to continue. Instead, the data grow quite rapidly toward 1, which is reached at 
Pt — 2 GeV, and possibly stay close to 1 at higher pr though experimental errors 
are too large to have a final say. This discrepancy is a mystery which cannot be 
explained up to now by current theoretical models, and needs further experimental 
investigation. 

The second and most interesting mystery appears when we look at negative 
rapidities, i.e., at the Au side. Fig. 2, compiled by A.Purwar [6], presents a summary 
plot of the pseudorapidity dependence of the Cronin effect for charged hadrons 
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Fig. 1. Top: i?cp measured by BRAHMS [4] at rapidity r\ e [0,3.2] for 2 centrality 
classes. Bottom: RdAu in minimum bias collisions, measured by BRAHMS [4] at 
rapidity rj G [0,3.2]. Note the striking difference of the forward rapidity RdAu '^'nd Rep- 
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Fig. 2. Summary plot of preliminary results from PHENIX and BRAHMS for integrated 
i?cp at 1 < pt < 3 GeV. Figure from Ref. [6]. 



integrated on 1 < pt < 3 GeV. As tlie rapidity decreases, the gold Bjorken's x 
increases. Then, the parton densities probed by the deuteron decrease and one 
would naively expect multiple scatterings and the Cronin effect to decrease, so 
that i?cp should tend to 1. On the contrary, experimental data show a steady 
increase as t] decreases, reaching R^^ ~ 1.8 at 77 = —2. This large value cannot be 
accounted for by standard nuclear effects like anti-shadowing [7] , or non perturbative 
string fragmentation which pulls hadron production slightly toward the Au side [8] . 
Moreover, the same trend at negative 77 appears also in J jil) production, ruling out 
a meson/baryon effect. 



2. pQCD and Glauber-Eikonal models 



The Glauber-Eikonal (GE) approach [9] to the Cronin effect treats multiple 2^2 
partonic collisions in coUinearly factorized pQCD, see Fig. 6, left panel. The low-pT 
spectra in p + A collisions are suppressed by unitarity. At moderate pr, the accu- 
mulation of transverse momentum leads to an enhancement of transverse spectra. 
At high j)T the binary scaled p -I- p spectrum is recovered: no high-p^ shadowing is 
predicted in this approach. 

The cross-section for production of a parton of flavour i, with transverse mo- 
mentum -pT and rapidity y ivL-p-V A collisions at fixed impact parameter &, is written 
as 

y b 



(PprdydPb 



i/plVi,PT 



(PpTdyi(Pb 



da'P 



d^prdyi 



(2) 



where where Ta (b) is the target nucleus thickness function. Hadron spectra are then 
obtained as a convolution with the appropriate fragmentation function Di^h{z, Q). 
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The first term in Eq. (2) accounts for multiple semihard scatterings of the 
parton i on the target nucleus; in the second term the nucleus partons are assumed 
to undergo a single scattering on the proton. The average parton flux from the 
proton, {xfi/p), and the parton-nucleon cross-section, da^^ , are defined as 



fi/p)vi,PT =^Yl i + Sij J ^y2a;i/i/p(a;i,Q2) ^ x2fj/N{x2,Q'^) j ^ 



•xdyi 
(3) 



where i is the Mandelstam variable and da/di are leading order parton-parton 
cross-sections in collinearly factorized pQCD. K is a, constant factor which takes 
into account ncxt-to- leading order corrections. To regularize the IR divergences of 
the single-scattering pQCD parton-nucleon cross-sections a small mass regulator po 
is introduced in the propagators, and Q = \/p^ + pfj2 is the scale of the hard 
process. Finally, we introduce a small intrinsic transverse momentum (fc^) = 0.52 
GeV^ to better describe the hadron spectra in p+p collisions at intermediate pt=1- 
5 GeV. The free parameters po and K in Eqs. (3) are fitted to hadron production 
data in p+p collisions at the energy and rapidity of interest. This allows to compute 
the spectra in p + A collision and the Cronin ratio with no extra freedom. For more 
details, see Ref. [9]. 

Nuclear effects are included in da , the average transverse momentum distri- 
bution of a proton parton who suffered at least one semihard scattering: 



da'^„ ... da'^ 



TA{b) X • • • X ^^Ta(6) e-'^'"(f«)^^('') . (4) 



X 



This equation resums all processes with n multiple 2^2 parton scatterings. The 
exponential factor in Eq. (4) represents the probability that the parton suffered no 
semihard scatterings after the n-th one, and explicitly unitarize the cross-section 
at the nuclear level. Unitarity introduces a suppression of parton yields compared 
to the binary scaled p + p case. This is best seen integrating Eq. (4) over the 
transverse momentum; da^^/dyd'^b ss 1 — e~'^ {po)TA(b) ^ i^r^ opacity x = 
a^^ {po)TA{b) <C 1, i.e., when the number of scatterings per parton is small, the 
binary scaling is recovered. However, at large opacity, X ^ !> the parton yield is 
suppressed: da^^/dyd'^b < 1 < a'^ {pQ)TA{b). This suppression is what we call 
"geometrical shadowing" , since it is driven purely by the geometry of the collision 
through the thickness function Ta- As the integrated yield is dominated by small 
momentum partons, geometrical shadowing is dominant at low pT- Beside the 
geometrical quark and gluon shadowing, which is automatically included in GE 
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Fig. 3. Cronin effect on pion production at Fermilab [10] and RHIC [11] at = 0. The 
solid curve is the GE computation. Theoretical errors due to the fit of po are shown as 
a shaded band around the solid curve. The rightmost panel shows the 0-20%/60-88% 
centrality classes ratio. 



models, at low enough x one expects genuine dynamical shadowing due to non- 
linear gluon interactions as described in, e.g.. Colour Glass Condensate (CGC) 
models, see Section 3. 

The GE model reproduces quite well both Fermilab and PHENIX data at 77 = 
(Fig. 3 left and middle). It also describe the increase of the Cronin effect with 
increasing centrality (Fig. 3, right). The agreement of GE calculations with y = 
data, and especially the centrality dependence, suggest that there is no dynamical 
shadowing nor CGC at RHIC midrapidity. 

To address the BRAHIVIS data at forward rapidity 77 « 3.2 [4], we would first 
need to fit pq and K in p + p collisions at the same pseudo-rapidity. Unfortunately 
the available pT-range pt ^5 4 GeV is not large enough for the fit to be done. 
Therefore, we use the parameters extracted at 77 = 0. The resulting Cronin ratio, 
shown by the solid line in Figure 2, overestimates the data at such low-px- However, 
the opacity xo = 0.95 might be underestimated, due to the use of the mid-rapidity 
parameters. To check this we tripled the opacity: the resulting dashed line touches 
the data at pt ^ 2 GeV, but predicts a strong Cronin peak at pt ?J 3 GeV. On 
the other hand, a similar strong peak would appear in the computation of Rep, in 
contrast with the data in Fig. 1. 



3. Colour Glass Condensate models 

The Colour Glass Condensate (CGC) is an effective theory for the nuclear gluon 
field at small- a: [12]. The valence quarks, treated as a collection of random colour 
sources p, radiate the gluon field A. At low-x the gluon occupation number is 
so large that the gluon field can be treated semi-classically and computed as a 
solution of Yang-Mills equation of motion in the presence of colour sources. The 
theory is characterized by a saturation scale Qs- Gluons with momenta lower than 
this scale are in the "saturation" regime: their density is so high that gluon-gluon 
fusion processes limit a further growth. At large momenta the gluon field is in the 
perturbative "parton gas" regime. A recently conjectured intermediate "geometric 
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d+Au^-h*+X 77 = 3,2 
\/s = 200 GeV Xa=0-95 
X BRAHMS h 




Fig. 4. Cronin ratio at rj — 3.2. 
The solid and dot-dashed curves 
are the GE model computation 
with X = Xo and x = 3xo, re- 
spectively. The dashed curve is 
a CGC computation by Jalilian- 
Marian [16] . Data are taken from 
[4]. 



[GeV] 



scaling window" may extend at momenta Qs < pt < Q^/Qo, where quantum effects 
from the saturation region further modify the evolution from the perturbative to 
the saturation region. 

Observables ©[^[p]] are computed as an average over the colour sources with a 
weight Wy[p] depending on the gluon rapidity y — log(l/x) as follows: (C'[^])y = 
/ Dp Wy[p]C'[^[p]]. Gluons at a given x are themselves colour sources for gluons at 
smaller x' < x. This evolution of the gluon field with x is captured by the so-called 
JIMWLK evolution equation [12]. 

If one approximates the proton as a dilute colour source, gluon production in 
p + A collisions can be explicitly written in a fcy-factorized form [13]: 



1 



(Pq^^dy l&TT^q\ J {2tt) 



r 1 e 



{uHo)u{x^)) 



Here, the proton's unintegrated wave function ipp is multiplied by the squared scat- 
tering amplitude on the target nucleus, expressed as a correlator of two Wilson 
lines. When the weight W is Gaussian, 



ri fy 

Wy[p] = A/'yCxp ^- / dy 



- y±) i 



>1 partonic 



gluon production in p + A collision can be interpreted as multiple 2- 
scatterings as illustrated in Fig. 6, center. 

To proceed further, one needs to choose a specific model for the colour sources. 
The McLerran-Venugopalan model assumes the gluon correlations to be local: A(ccj^— 
V±) = Qsix)6{xj_ — yj_) where Qs{x) = a^^^xG{x,Q^). Quantum evolution of 
the gluon field is included naively in the x-dependence of Qs only. An estimate from 
Ref. [13] of the Cronin effect on gluon production is presented in Fig. 5, left. It 
qualitatively agrees with experimental data at y = 0. However, at j/ = 3 it predicts 
a rightward shift and an increase of the Cronin peak (analogously to the GE model). 
It is disfavored by the Rep data, though not incompatible with RdAu- 

On the other hand, at larger rapidities (i.e., at small x) the system undergoes 
a long quantum evolution, and enters the saturation regime. Thus, we should solve 
the JIMWLK equation. In the Gaussian approximation, a self-consistent solution 
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Fig. 5. Cronin effect in the MV model with naive quantum evolution (left) and with 
full evolution in the "deep saturation" limit (right). Figure taken from Ref. [13]. 

in the saturation regime exists, but with non-local gluon correlations [14]. In other 
words, the gluon density has become so high that gluons are no more bound to 
single nucleons, but correlate over macroscopic distances. The resulting Cronin 
ratio is highly suppressed over the whole pT range, in qualitative agreement with 
Rcp{y = 3). However, it disagrees both with the forward RdAuiv = 3) and with 
mid-rapidity data. 

A word of caution: at large rj valence quark production is important, if not dom- 
inant, as demonstrated by the excess of positively charged hadrons over negatively 
charged ones in p + p collisions at ry = 3.2 anf pr > 2 GeV [15]. Therefore, these 
computations should be considered only as illustrative of the physics included in 
the CGC approach. However, even the more quantitative computation of Ref. [16], 
who takes into account valence quarks and their fragmentation into charged hadrons, 
fails to reproduce the pt > 1.5 GeV of RdAu at y = 3, see Fig. 4^. 

4. Where does pQCD meet the CGC? 

At mid-rapidity, GE and CGC models give similar results for the Cronin effect. 
Does this mean that to some extent they are equivalent, even if GE models de- 
scribe multiple 2^2 scatterings and CGC multiple 2^1 scatterings? The standard 
answer is yes. Consider the single scattering terms depicted in Fig. 6, right. The 
single-inclusive cross section in pQCD is obtained by integration over the unob- 
served parton rapidity y' . In the limit of large px ^ Qs and asymptotic energy 
PT/\/se~y <C 1, one obtains in the MV model dapQCD ~ dupQCD oc l/py [18]. 
This is usually read as a proof that the single-inclusive cross section in the MV 
model reduces to pQCD in the high-pT limit. Hence the GE and the MV models 
are equivalent. 

However, the above equivalence can hold only in a restricted kinematic domain. 

This is clearly seen if we compute the average Bjorken's x's probed in p+p collisions. 
In the 2^2 kinematics, they depend on the rapidity of both the observed and the 

-"^A very recent computation [17] was more successful. It is based on a parametrization of 
the onset and rapidity dependence of quantum evolution which mimics CGC effects, and models 
valence quark scatterings in a different way than in Ref. [16]. 
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GE: multiple 2^2 



CGC: multiple 2^>1 



dy- y 




y =r^r y 

PT>Qsat — 



Au ^Au 

Fig. 6. Left and middle: Multiple scattering processes included in the GE and CGC 
models, respectively. Right: Illustration of the equivalence of pQCD and the MV model 
at high px discussed in Ref. [18] . 



unobserved parton: 



.±y' 



In 2-^1 parton processes, the observed particle fixes them completely: 

.2^1 _ PL(.±y 



-1.2 



As shown in Fig. 7, the 2^1 and 2^2 processes probe the target at quite different 
Bjorken's x\ The reason of this difference is that in the 2— >2 processes one has to 
integrate over the unobserved particle, over a quite large phase space Ay' « 10. 
When y is large, most of these particles are produced at small or negative rapidity, 
which compensate the decrease in {X2) expected when the rapidity of the observed 
particle increases. 

Given Fig. 7, one is tempted to ask "Which of the two models is wrong: GE or 
CGC?" In my opinion, this is the wrong question. The existence of 2— >2 hard scat- 
tering processes at mid-rapidity has been experimentally proved by the observation 
of back-to-back hadrons at large pr in p+p and p + A collisions [21]. On the other 
hand, there is no reason to doubt about the existence of 2-^1 processes. So a more 
pertinent question is "What is the interplay between 2-^2 and 2— >1 processes? Is 
there a region where one is dominant over the other?" 

We may have a clue at the answer from the following zeroth-order argument. 
Neglecting the quarks for ease of notation, the single inclusive cross-section for 2^2 
processes may evaluated as: 



= I dy xiGixi ) X2G[X2) 

dy J ' dt 



da 



« Ay' X xG{{xi)2-,2) XG{{X2)2^2) 

where Ay' is the width of the integration interval over y', and (a:i^2)2-+2 can be read 
off Fig. 7. In the 2^1 case we have: 



daPP^sX 
dy 



xG{{xi)2^i) xG{(x2)2^i) 



da 
di 



where the value of (a;i.2)2^i is shown in Fig. 7. At mid-rapidity, (2:2)2^1 ~ (2^2)2^2, 
then 2-^2 processes are dominant because of a larger phase-space (Ay' w 10). At 
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Fig. 7. Average (xi) and 
(12) in p+p collisions at RHIC 
energy, for 2— >1 processes 
(left, taken from Ref. [20]) and 
2— »2 processes (right). 
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very forward rapidity, (a;2)2-^i ^ {x2)2^2\ hence, 2— s-1 processes are dominant 
because of the small-x growth of parton distributions (the same argument can be 
used to analyze backward rapidities). 

The transition between these two regimes is poorly known even in p + p colli- 
sions, and should be more carefully studied, both theoretically and experimentally. 
This requires a formalism which is able to describe both processes in a unified 
framework. CoUinear factorized pQCD is not suited, as kinematics forbids 2^1 
processes if incoming partons have zero transverse momentum. On the other hand, 
kT factorization can describe 2— >1 processes, and the possible production of a sec- 
ond or more hard particles is hidden in the unintegrated distribution functions for 
the incoming partons [19,20]. It is thus very suited for the problem at hand. 

5. Conclusions 

The pQCD based GE model is able to quantitatively describe the Cronin effect 
at midrapidty over a broad ^/s — 20-200 GeV energy range, indicating absence of 
CGC. However, at large rapidity it predicts an increasing Cronin effect, which is 
unfavoured by experimental data. Only incoherent 2-^2 multiple parton scatterings 
are included in the GE model; adding coherence effects would suppress the Cronin 
peak, but still predict RdAu ^ Rep 1 at large pr [22]. On the other hand, CGC 
models describe 2^1 multiple coherent parton scatterings, and suggest at RHIC 
a transition from a non saturated region at mid-rapidity, to a saturated region 
at forward rapidity. Here, nonlocal gluon correlations created during quantum 
evolution suppress below 1 the otherwise increasing Cronin ratio. This predicted 
plateau at high-pr is the hallmark of CGC. Unfortunately, the discrepancy between 
RdAu and Rep BRAHMS data, which is not present in any of the discussed models, 
does not allow to conclude whether this picture is correct or not. 

At mid-rapidity GE and CGC models are usually thought to be equivalent, 
even though they describe different partonic subprocesses, and probe quite different 
Bjorken x's at forward and backward rapidity. Simple kinematic considerations 
suggest a dominance of 2^2 processes at mid-rapidity, and of 2^1 processes at 
forward and backward rapidity. The interplay and transition between the two is 
poorly known, and needs to be studied more accurately already at the p + p level, 
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then extended to p + ^ collisions. The appropriate tools are already on the market: 
fcT-factorized pQCD on the theory side, and the STAR and PHENIX detectors 
on the experimental side. Azimuthal two-hadron correlations and their rapidity 
dependence arc an especially suited observable for this purpose. 

The steep increase of Rep at negative rapidity adds a new dimension to the 
problem. Indeed both GE and CGC models expect the Cronin effect to be more 
and more reduced at 77 < 0. This is due to the decrease of the nucleus opacity XA{y) 
in one case, and of the saturation scale Qs{y) oc in the other. Standard nuclear 
effects like anti-shadowing are too weak to explain the magnitude of the negative 
rapidity Cronin effect, which requires a new theoretical explanation: either a ne- 
glected effect in the GE and CGC models, or a new piece of physics. 

Acknowledgements. I would like to thank M.Gyulassy, P.Levai, A.Purwar, 
J.W.Qiu, A.Szczurek and X.N.Wang for inspiring discussions. 
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